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7. Introduction 


THE study of the transmission of light through a cloud of particles is of 
practical importance, since it is connected with natural phenomena such as 
the opacity of clouds and smokes, and the visibility in fog and rain. So 
far, such a study has been restricted only to those cases where the particles 
are very small. Lord Rayleigh (1899), for instance, has used his theory of 
the scattering of light by small particles to determine the attenuation of the 
transmitted light. There is, however, no satisfactory simple theory, either 
for transparent or opaque particles of bigger size, distributed in a transparent 
medium. Mallock (1919) has attempted to give a theory of the transmission 
of light through a cloud of water droplets, but his theory is defective in 
many respects. He has assumed the droplets to be opaque, which is not 
true. Also, he has treated the problem on a geometrical basis, omitting the 
effect of diffraction, which is not justified even with opaque particles. 
Thirdly, his application of the theory of probability to the geometrical method 
is also not correct.* 


In this paper, a general theory is developed of the propagation of light 
through a cloud of spherical particles, randomly distributed in a transparent 
medium. Both the cases when the particles are opaque, as well as when they 
are transparent are considered from the point of view of the wave-theory 
of optics. It is shown that in general, there is always an attenuation of the 
transmitted beam, which is of the exponential form. The numerical value 
of the attenuation coefficient has also been calculated in a number of cases. 
In particular, Mallock’s formule are corrected during the course of the 
development of the theory. 


2. Geometrical Treatment of the Problem 


In this section, we shall develop a theory of the effect on the trans- 
mitted light of a number of opaque spherical particles distributed at random 





* Richardson (1919) has also given a similar theory, but he too has taken the droplets to be 
Opaque, and has not taken diffraction into account. 
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in a transparent medium, purely on the basis of geometric ideas. The case 
in study may be illustrated by the transmission of light through heavy smoke, 
where the particles composing the smoke can be considered to be opaque. 


Suppose that a beam of cross-sectional area A traverses a column of 
the smoke of length /, and that N is the number of particles per unit volume 
of the medium. Assume further that the particles are all of the same size, 
and have a radius a. On a geometric basis, the effect of each particle would 
be to cut out a portion of the transmitted beam of area 7a*, and thus reduce 
the energy content by an amount corresponding to this area. In this way, 
each particle reduces the energy content of the beam by a certain fraction, 
and it would thus appear that if the number of particles my be such that 
Npra? = A, then the beam would be completely cut out. However, it is not 
so, on account of the fact that some of the particles would screen those behind 
them, and thus increase the chance of a portion of the direct radiation 
coming through. It is therefore a question of probabilities to determine 
what fraction of the energy is transmitted by the particles in the medium. 
In our case, the number of particles, 2, is evidently = NAI. 


The problem at hand is identical with one in which n disks, each of 
radius a, are thrown at random on an area A, and it is required to find the 
probable area covered by the n disks. Let P, be the probable fraction of 
the area covered after m disks are thrown, and let p be the fraction of the 
area covered by a single disk when alone, i.e., p =7a?/A. Now, on adding 
one more disk, it may fall on the area already covered, or on the area un- 
covered. The probability of its falling on the empty area is evidently 
(1 —P,,). Hence, the probable fraction of the area covered after (m+ 1) 
disks are thrown is 


Pri = Pmt p(l — Py) =P (1 —p) +p (1) 
Now, P, =p, so that 
P, =p(1—p)+P, 
P; =p (1—p)*?+p(l1—p)+ Pp, etc., giving 
P= 2 1 —py*+ 6 — pr +--- + 
=p [i—(1—p)*)/[l-— (— p)] = [1-(—p)". (2) 


Thus, the probable fraction of the area not covered after w disks are 
thrown is 


Q= 1— P= (1— p)”. (3) 
Coming back to the problem of the transmission of the beam of 
light, the probable fraction of its area A, not cut out after it has encountered 
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n particles is given by (3). Hence, the intensity of the transmitted beam is 
I,= I; (1 — na*/A)A™, (4) 
where I;= intensity of the incident beam. Since za*/A is a small quantity, 
we get 
Wh. (5) 
Thus, it is seen that the intensity of the transmitted beam is never 
zero, but that it diminishes exponentially with increase of the thickness of 
the medium. Also, from the relation Q,=—(1—p)*, it is seen that the 
effects of the particles are multiplicative, each one diminishing the intensity 
in the ratio (1—p):1. This result is important, since we have started 
with the assumption that the effect of each particle is subtractive, i.e., each 
subtracts a portion of the energy from the incident beam, corresponding 
to its area. However, on account of the random distribution of the parti- 
cles, the effect becomes multiplicative. 


Incidentally, it may be noted that if n= 1/p, i.e., if the total area of all 
the particles is equal to the area of the aperture, then Q, = (1 — 1/n)*= e". 
Hence, the probable fraction of the energy coming through is 1/e. Mallock 
has wrongly assumed the value 1/2 for this quantity without any proof. 


The actual value of the attenuation coefficient derived in this section 
is not correct, for, as is shown in Section 3, diffraction effects produce a doubl- 
ing of its value. But, before proceeding to consider these, it will be inter- 
esting to discuss what the attenuation would be, according to geometric 
ideas, if the particles were not all of the same size. In this case, we must 
expect a distribution of size among the particles. As a general case, suppose 
that the number of particles per c.c., having a radius between a and 
a+ da is 

dN = f (a) da, (6) 


where f(a) is the function giving the law of distribution. 


Now, since all the particles are not of the same size, the application 
of the theory of probability has to be modified. Let p,, pe, p3,°** be the 
fraction of the area covered by the first, second, third, etc., particles, when 
alone. These can be equal or different. Then, we have the general relation 
similar to (1), viz., 


Preit ae P,, + Pm+i (1 Pas ®) an P., (1 aes Pmvi) + Pm+1 (7) 
Hence, 


Pi= pr, Q, = (1— p,) 
P.= p, (1 — pz) + pro, Q. = (1 — p;) (1 — pa) 
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Ps = P; (1 — pz) (1 — Ps) + Po (1 —Ps)+ Ps, Qs=(1— py) (1 — pe) (1 — py), 
etc. 


Thus, every one of the quantities Q, is symmetrical in p,, p2, --- so that 
the order of the quantities p,, p2,"* +P, is immaterial. Hence, we may regard 
each particle as producing an effect independent of the rest, and reducing 
the energy of the transmitted beam in the ratio (l1—p,):1. Hence, 
expression (4) can be generalised, and written as 


-Ep ” 
1/1; = (1 — py) (1 — pe) + "(1 — py) = € 
If we denote by a, the radius of the nth particle, p,,.=a,?/A. Also, on 


account of the distribution (6), the summation 2p, can be written as 
1 


Sna,2/A = ml fatdN. (9) 
Hence, 1,/I, = exp [— al fa*dN] 
= exp [— zla®N] (10) 


where a? is the mean square of the radius. 


As a particular case, we may take the Maxwellian law of distribution 
given by f(a) = a*e-“, or 
dN = a*e™ da. (11) 


In this case, a? = Jn 


a*, where @ is the mean radius, so that 
I,/I; = exp (— 3n?4?NI/8). (12) 


This expression can also be put in terms of the most probable value 
of the radius, a,,, as 


i 


1,/I,= exp [— 3x a? NI/2]. (13) 

So far, we have been dealing with spherical particles. This assump- 
tion is, however, unnecessary, for the shape of the obstacle is unimportant. 
The really important quantity is the area of the beam that it cuts out. If 
we represent the area of a particle by a, then, p = a/A, and (5) can be written 


in the form 
1,/]; = exp [- alN]. (14) 


In the case, where the size of the particles is not the same, 
1,/I;= exp [—/ fadN]. (15) 
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3. Transmission according to the Wave Theory 


It is well known that a small opaque circular obstacle placed in the 
path of a beam of light produces a diffraction pattern, whereby part of the 
energy is diffracted in directions away from the forward one. On account 
of this, it is evident that the reduction of energy in the forward direction 
will be greater than that due to the mere cutting off of a portion of the wave- 
front by the obstacle. In the following discussion, it will be shown that the 
reduction is exactly double of what would be expected on the geometrical 
theory. 


In this connection, it is important to consider what we mean by the 
intensity of the transmitted beam. The aperture of the incident beam is 
always finite, and will have a definite area A. Consequently, this will 
produce a diffraction pattern of its.own. If we take the amplitude of the 
incident wave as unity, the intensity in the exact forward direction would 
be A2/A2, and the pattern would extend over a solid angle of the order of 
\?/A. We may therefore speak of two quantities as representing the inten- 
sity of the transmitted beam: (a) the intensity in the exact forward direction, 
and (b) the total intensity in the diffraction pattern, which is in general 
comprised within a small solid angle. If the beam is brought to a focus on 
a photographic plate and we take a microphotometer record of the intensity, 
then the height of the peak would represent the quantity (a). If, on the 
other hand, we focus the beam on a photocell, the current produced would 
represent the quantity (6). The distinction between the two is necessary, 
since the reduction due to the existence of an obstacle is not the same for 
both, as is shown below. 


We will now consider the general method of finding the reduction in 
the forward intensity. Let U be the amplitude of the diffraction pattern 
produced by the aperture A. Then I = /| U|*dQ, where @ is the solid angle, 
represents the total energy in the diffraction pattern, the integration being 
done over all the solid angles where |U |? is finite. On the other hand, 
I,= | U,|? represents the intensity in the forward direction per unit solid 
angle, if U, is the amplitude in this direction. 


Now, suppose that an obstacle of area a be introduced in the path of 
the beam. Then,.this would produce a diffraction pattern of its own. Let 
u be the amplitude in the diffraction pattern due to an aperture comple- 
mentary to the obstacle. Then, the amplitude in the pattern due to the 
obstacle is, by Babinet’s principle equal to —u. The two patterns due 
respectively to the aperture A, and the obstacle a will be superposed, the 
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second of which would obviously cover a larger solid angle. In this case, 
the total energy in the transmitted beam is 

l’=/|U —ul* dQ, (16) 
where the integration is to be carried out over the same range of solid 
angles, as was done for I. On the other hand, the intensity in the exact 
forward direction is 

Ig’ = |Uo— uo’. (17) 
From these formule, the reduction in intensity either of I or of I, can be 
calculated. 


Coming now to the particular case in which we are interested, assume 
for simplicity that the obstacle a is spherical, and has a radius a. The 
diffraction would now consist of two parts, (1) the pattern due to the 
aperture A, covering a solid angle of the order of A?/A, and (2) the pattern 
due to the obstacle, covering a much larger solid angle of the order of A*/a. 

Evidently, I, = A?/A*, and I,’= (A — a)?/A*, so that the ratio I,'/I,= 
(1 — 2a/A), since a is small. Thus, the intensity in the exact forward 
direction is reduced in the ratio (1 — 27a?/A): 1, or double the reduction 
as given by the geometric theory. 


In order to evaluate I’, assume that the pattern due to aperture extends 
over a solid angle 2,; 2, is ~ A?/A. Then, 
I= f|U|?dQ =A. (18) 
25 
Take now the whole pattern due to the aperture and the obstacle, and 
assume it to extend over the solid angle 2. Then the total energy in this is 


f\U —ul* da. (19) 
Q 
But, outside 25, U =0, so that the above expression can be written as 
f\U —ul?dQ+ ful? dQ. (20) 
Q, Q —2, 
Hence, 
V=/|U —ul*dQ =/\U — u|?dQ —f |u|? dQ. 
2, Q Q — 25 


Now, 2, is € 2, so that the energy of the pattern due to the obstacle con- 

tained within Q is very negligible. Hence, f |u|? dQ= f |u|? dQ, which 
2-2 Q 

must be equal to a, by Babinet’s principle. Also, /|U — u|* dQ will be equal 
Q 


to (A —a), so that 


Y=A — 2a. 
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Thus, l/l = (1 — 2a/A), (21) 


i.e., the total intensity is also reduced in the same ratio as the forward 
intensity, the reduction being double that given by the geometric theory. 


A physical interpretation can also be given for this double loss in the 
total transmitted intensity. The cause has already been explained as due 
to the effect of diffraction. Now, by Babinet’s principle, the energy in the 
diffraction pattern due to an opaque disk is the same as that in the pattern 
due to an aperture of the same size, and hence equal to a. In this way, the 
forward wave loses an amount of energy a due to diffraction, in addition 
to an equal amount lost due to the blocking by the opaque obstacle. This 
is exactly what is shown by the formule above. 


However, it must be noted that we have assumed in this derivation 
that the energy of the diffraction pattern due to a, contained in the solid 
angle embraced by the pattern of the transmitted light is negligible. This 
will be true, only if a < A. However, if this is not the case, then the 
additional loss of energy will not be a, but less. In general, therefore 


I'/I =(1 — ca/A), (22) 


where c = 2, only if a is small compared to A. If the obstacle is of dimen- 
sions comparable to those of A, then the pattern due to it will be quite small, 
and we have to take c nearly equal to 1, as is found from geometrical consi- 
derations. This is quite natural, since geometrical optics applies rigorously 
when large sizes are taken into consideration. However, in the exact forward 
direction, the reduction in intensity is always double that given by the 
geometric theory. The double loss of energy, and the conditions determin- 
ing its occurrence do not seem to have been pointed by anybody so far. 


We will now consider what the effect of a cloud of particles would be 
according to wave theory. Assuming that the incident wave is represented 
by u= exp (2zict/X) at x=0, which corresponds to the beginning of the 
cloud, the wave at a distance x within the medium can be represented by 


2ni 
—7* (et — wx) 
u—-e**e A - 


where p» is the refractive index, and k, the absorption coefficient. Taking 
now a thin layer of the medium of thickness d/, the number of particles 
in this thickness per unit area of cross-section will be Nd/. As already seen, 
the effect of each particle would be to reduce the amplitude by a. Hence, 
the wave, after passing the thickness d/ would be represented by 


2ni 
a~Meapes oO: 
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Comparing this with the above equation, we get »=1, and exp (— kdl) 
= (1— Nadl)= exp (— Nadl). Hence, the absorption coefficient for the 
amplitude is k= Na, and that for intensity will be double this, namely 
2Na. Hence, the intensity of the beam after passing through a length / of 
the medium is 

p= 1, ene! (23) 
where I; is the incident intensity. For spherical particles of radius a, this 
gives 

L/T, = e72e@n?, (24) 
This expression differs from that given by the geometric theory only in 
that there is an additional factor 2 in the exponent, the reason for which 
has already been explained. 


4. Passage of Light through Transparent Spheres 


This case is of practical importance in connection with the study of 
the transmission of light through a cloud of water droplets, as in fog, 
cloud or rain. Mallock has wrongly supposed tnat one can take the water 
particles as opaque. This is, prima facie, not true, since the droplets are 
by their very nature transparent, and not opaque. Consequently, portions 
of the wave-front transmitted through the droplets must also be taken 
into consideration. The portions passing through the droplet will, however, 
be retarded with respect to the rest of the wave-front, the extent of the 
retardation depending on the distance each portion has passed through the 
drop. Changes in amplitude will also occur in the wave-front on account 
of the convergence of the rays passing through the drop. These features 
must be taken into account in a proper theory of the phenomena. 


The actually observed facts also show that the drops cannot be regard- 
ed as opaque. The work of Barus (1907, 1908), and of others has shown 
that the transmitted light is often highly coloured and that this colour is 
complementary to the colour of the light scattered by the particles them- 
selves in the forward direction. Also, these colours are found to undergo 
periodic changes with increase in particle size, reappearing for values of the 
radii in the ratio of natural numbers. These phenomena indicate that the 
light transmitted through the drop must be taken into account. 


The most rigorous method of attacking the problem would be to use 
the electromagnetic theory of light, and to determine the effect of a 
dielectric sphere on the field. Such a theory has been given by Mie (1908), 
Debye (1909), Rayleigh (1910) and Bromwich (1920). However, this method 
does not give any simple expression for the intensity of the diffracted and 
the transmitted light, if the size of the droplets is not small compared with 
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the wave-length of light. The nature of the phenomena could only be 
studied by numerical computation, which again is prohibitively tedious. 
It is therefore desirable to give a theory based on optical principles, which 
would yield expressions that can readily be applied. 


We shall now attempt to give such a theory for the diffraction of light 
by a transparent sphere. The theory will be developed to include also the 
diffraction of light in directions inclined to the forward one. As already 
said, one must take into account both the amplitude and phase changes 
that occur in the transmitted wave-front. However, if we assume that the 
difference in refractive index between the sphere and the surrounding medium 
is small, then we can as a good approximation, take that the effects of 
refraction are negligible. In this case, we can neglect the changes in 
amplitude, and take only the phase-changes into consideration, the path 
retardation of any portion of the wave-front being calculated on the 
assumption that it has passed straight through the drop. 


T. A. S. Balakrishnan (1941) has developed a theory for the diffraction 
by a single droplet based on similar ideas. In this paper, the theory is 
worked out more fully, and it is extended to the consideration of the 
transmission and forward scattering by a cloud of transparent spherical 


droplets. 
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Fic. 1. Passage of Light through a Sphere 


5. Derivation of the Formule 


In this section, the formule for the diffraction of light by a spherical 
droplet is derived, for the case when the amplitude changes can be neglected, 
and the wave-front can be supposed to have passed straight through the 
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drop. Let P,RP,’ be a plane wave-front which, after passing through 
a droplet of centre O, and radius a, emerges out as P,R’P,’ (Fig. 1). The 
phase of the emergent wave-front will evidently not be uniform throughout, 
but will lag behind in the portions which have passed through the drop. 
Let RR.R’ be any ray, and let the angle between OR, and OO’ be @. 
Then, this ray suffers a retardation of 2a (u— 1) cos @ with respect to the 
portions of the wave-front passing outside the drop. Thus, the plane 
wave-front, on the whole, would be distorted into one having a concave 
part in the centre. 





Our idea is to find the Fraunhofer diffraction pattern of this emergent 
wave-front, i.e., we wish to find the intensity in a direction making an 
angle 4 with the incident direction. For this, we proceed in a manner 
analogous to the usual theory of diffraction by a circular aperture (vide 
Max Born, Optik, pp. 155 et seq.). ‘aking polar co-ordinates, if (r, a) be 
the co-ordinates of an element at R’ with respect to O’ as origin, and the 
horizontal radius as the initial position of the radius vector, then the 
diffracted amplitude in the direction ¢ due to this element is 


sin 5 (ct—2(u— 1) acos 9+ rcosa sin ¢} r dr da, 


where the amplitude of the incident wave-front is assumed to be unity, 
and A is the wave-length of the light. Substituting a sin @ for r, and 
integrating between the proper limits, we get that the amplitude in the direc- 
tion @ due to the portion of the wave-front which has passed through the 
drop is 





n[2 a 
a* a: ae 5 . ; 
X, ae | sin + {ct— 2 (u— 1) acos 6+ asin @ cos a sin 4} 
0 0 sin 6 cos @ dé da. 
Putting x= 2mct/A, € = 4m (u— 1) a/A, and n= 27 sin ¢ a/A, 
n[2 2n 


no 





x, « ‘ sin (X — € cos 0+ 7 sin @ cos a) sin 0 cos 0 dO da. 
5 
Integrating with respect to a, this can be written as 
nf2 
27 a? . , . 
x, = 5 f Jy (7 sin 8) sin (X — £ cos 8) sin 8 cos 0 d8. (25) 
if 


The effect of the portion of the original undisturbed wave-front covered 
by the drop is obtained by substituting » =1 in the above expression, 
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which makes = 0, and gives the value 
2 


n/ 
a 
X, = —— f Jo (m sin 4) sin 6 cos 6 dé 
0 


= aie hin - sinX, (26) 
The amplitude X of the diffracted wave due to the complete wave-front 
which has traversed the drop is the sum of the amplitude due to the 
portion which has passed through the drop, and that due to the portion 
which has passed outside it. The former is given by X,, and the latter by 
subtracting X, from the effect due to the undisturbed wave-front, which we 
may call X,. In the region outside that covered by the source, X, is 
evidently zero, so that 
X= X,— Xz. (27) 
In the forward direction, i.e., for 6=0, Xp is finite, so that the resultant 
amplitude, X’, in the forward direction, is 


X’ = X_+ X,’— Xe’, (28) 
where X,’ and X,’ are the values of X, and X,, when ¢=0. For this 


particular case, since ¢=0, y vanishes, and the integral in the expression 
for X, can be easily integrated, giving 





n{2 
X,’ = ie f sin (X — £ cos 9) (€ cos 4) d(€ cos 4) 
0 
és ei sin X [€ sin + cos €— 1] + se cos x [F cos:¢'— stag}. 
_ 7a* i) 


Also, X,’ = sinX, since Lt 


A 7-0 
X’ = X, + 2na2/A-sin X -(cos €/€?+ sin €/€ — 1/2— 1/€*) 
+. 2na?/A-cos X (cos €/€ — sin €/€?). 
Now, if A be the aperture of the beam, then X, is clearly equal to 
A/X. sin X, so that 


, ‘ , 
x= Ast Peat {oone 4 ES - zt] 


= 4, so that 





RUF — 2 & 
A cosX [27a? (cos é _ sin € ] 
ee Ar etl (29) 
Now, the quantity 27a?/A is small, (since A can be made as large as we 
please), so that its square can be neglected. Hence, the intensity in the 
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forward direction due to the wave-front that has passed through a droplet is 


ee — Sfr— + Soret ey es Se rey ee 


~ 2 A l2°8 & g 
This expression gives the reduction in intensity due to a single droplet. In 


order to find the effect of a cloud of droplets, we proceed in a manner 


similar to what we adopted for a cloud of opaque particles. Let the incident 
ey : cet , 
wave be represented by u = sin I< (ct — nx), the origin being at the 


boundary of the cloud. Then, from (29), the amplitude of the wave after 
passing through a layer of the cloud of thickness dl is 


2 
u’ == [1 — Ndl 2na® f()] sin + (ct — dl) — Ndl 2za’g (é) cos = (ct —dl), (30) 
where g (&) is the expression (sin €/€*— cos €/é). From the wave-theory, 
we know that wu’ must be of the form 
u’ = e* sin - (ct — pdl) (31) 


Putting the expression (30) above in this form 
u’ = [{1 — 2xNa? f (£) dl}* + {2nNa? g (€) dl}*}"? sin {= (ct— dl) —« \, 
where tan ¢€ = 27Na’g (€) dl/{1 — 2aNa? f (€) dl}. 


Since both the quantities 27Na?f(é) dl and 27Na?g (é) dl are small com- 
pared to unity, this may be written as 


y= [1 — 2aNa* f(é) dl] sin (2m (ct — dl) — ¢} 
== @- 2HNG*/(g)4? cin {z (ct — dl)—el, 


where tan ¢ = 27Na? g(é)dl = «. 
Comparing this with (31), we see that k = 27Na?f(é) and » = 1 + Nda’g (é). 
Since the absorption coefficient for intensity is double that for amplitude, 
we get that the intensity after passing through a length / of the cloud is 

I, = I, exp [— 4za°NI f(€)], (32) 


where I, is the incident intensity. 
6. Colours of the Transmitted Beam 


Equation (32) above gives the formula for the attenuation of the directly 
transmitted beam. Ordinarily, for particles whose radii are larger than 
the wavelength of light, é is large, so that terms in 1/&* can be neglected. 
(32) then goes over into the form 


1,/Ig = exp [— 4xa?N/ (1/2 — sin €/8)]. (33) 
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This expression clearly shows that the attenuation of the transmitted light 
is different for different wavelengths, the attenuation coefficient depending 
on é, and hence on A. Taking particles of a fixed size, the value of € will 
evidently depend on the wavelength A of the light, and from (33) it is clear 
that the intensity will be a maximum for those for which sin € is + 1, and 
a minimum for which sin € is —1. Now, the visible spectrum comprises 
the range of values 0-42 to 0:7 of A, so that if the drops are small enough, 
there will be only one maximum, and the colour of the source will be due 
to the predominant effect of the range of wave-lengths in the neighbourhood 
of this maximum. This explains the reason why the transmitted light is 
vividly coloured. In Fig. 2, f(€) is plotted against ¢, and it shows the manner 
in which f (€) undergoes maxima and minima. 


°, I | \ tin 
| 





a : 20 
(4) 40 

& 

a J 








Fic. 2. Graph of f (¢) against ¢ 


As the size of the droplet is changed, the value of A for which the 
maximum occurs will vary, and consequently, the colour of the transmitted 
beam will also vary. Since the varying function is sin é, we should expect 
a cyclic change of colours as é increases; the colours will repeat themselves 
for values of € which are in the ratio of natural numbers, i.e., for values of 
the radius a which are in this ratio. This is in accordance with Barus’ 
observations (1907). 


If the size of the droplets is further increased, the range of values of é 
may comprise two or more maxima, and also the difference between the 
maxima and minima of intensity will not be pronounced, since (sin £)/£ 
becomes intrinsically small. Hence the colours will become less and less 
vivid. In the limiting case of very large values of £, sin ¢/£ becomes small 
compared to 4, and (33) reduces to 


I/Ip = exp [— 2na*N]]. (34) 
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There is no term depending on A, so that no colours will be exhibited. It 
is also interesting to notice that this is identical with the intensity in the 
forward direction given by a large opaque sphere of the same size. The 
following simple physical argument shows that this must be the case. When 
the particles are large, the path retardation between the central ray, and 
the rays passing outside the drops will comprise a number of wavelengths. 
We can therefore divide the portion of the wave-front which has passed 
through the droplets into half-period zones, of which the outer ones will 
practically cancel out. The resultant effect will only be due to the innermost 
rings, the area of which relative to the total area of the wave-front that has 
passed through the drop is smaller, the larger the drop. Thus, the forward 
intensity must approach the value for an opaque disc. In fact, the vanish- 
ing of the term (sin £)/€ mathematically represents this fact. 


7. The Intensity of Forward Scattering 


As already mentioned, we can employ our theory to consider the nature 
of the light scattered in the forward direction by the droplets. This would 
comprise the region just outside the region covered by the source. Here, 
¢ is small, but the effect of the incident wave is zero. The amplitude is 
therefore X = X,’— X,’. Substituting for X’, and X’., 


x = 200 sinx [2958 + SBF 5 2] + 2a cosx [Ef SF]. 5 


A 2 ¢ A é Ed 
As before, neglecting terms in (1/&*) as an approximation, we get the intensity 
as equal to 
(2ma*/A)? (1/4 — sin €/£). (36) 

This expression refers to the intensity diffracted by a single droplet. But 
the light reaching it, and the light diffracted have both to travel through the 
fog, and the total path will not differ from /, since ¢ is small. Hence, the 
expression (36) has to be multiplied by the expression 

exp [ — 4na?NI/ (1/2 — sin €/8)]. 
Also, since we have assumed a random distribution of the droplets, there 
will be no phase correlation between the light scattered by different particles. 
Hence, the light diffracted out of the cloud in the forward direction, per 
c.c., is I; = N (2ma?/A)*- (1/4 — sin €/£). 

exp [— 4za?N/ (1/2 — sin €/€)]. (37) 
This is of the form Bxe-, where x = N (1/4 — sin €/£). The correspond- 
ing expression for the transmitted light is I,=Ije-*. Now, the former 
expression is zero when x = 0, and reaches a maximum for x = 1/C, after 
which it diminishes. The latter is a maximum at x =0, and steadily 
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diminishes as x increases. Hence, if the maximum value of x, upon which 
the colours depend, is less then 1/C, then corresponding to it I¢ will be a 
maximum, while I; is a minimum. Similarly, for the minimum value of x, 
the reverse will be the case. Thus, the colour in the region immediately 
surrounding the source will be complementary to that of the source itself. 


However, it is to be noted that for this phenomenon to be prominent, 
the maximum value of x must be less than 1/C, or N/ must be less than a 
quantity which depends on the size of the droplets. In other words, if the 
density of the fog is fixed, there is a limit to the depth of the fog beyond 
which the colours are not complementary, but will tend to be more and more 
identical. The phenomenon of complementary colours can best be observed 
with fairly thin columns of cloud, as for instance when the sun is seen 
through the puffs of steam emitted by a locomotive. 


In the limiting case of very large drops, the expression (36) for the 
forward intensity reduces to (7a*/A)*. This is again identical with that for 
an opaque obstacle of the same size. 


8. Concluding Remarks 


The theory developed in this paper for the transmission of light through a 
medium containing a cloud of particles cannot be expected to be strictly valid 
for very small sizes of the particles, since Huygens’ principle, and the idea 
of complete wave-fronts, on which it is based cannot be applied when very 
small dimensions are concerned. The correct procedure in these cases is 
to apply the theory of scattering of light, as for example, the theory developed 
by Mie and others (/oc. cit.). These however do not yield any simple 
expressions ; but for very small droplets, the theory of Rayleigh can be used. 
In this case, the attenuation has actually been worked out by Lord Rayleigh 
(loc. cit.) and can be represented by 


12875 = (u?—1)? af 

3 (uw? + 2)? x: 
In the case when (u— 1) is small, (u? +2)? can be put as equal to 9, 
so that 





1,/Ip= exp [ _—NI (38) 


I)/Ig= exp [— 12875 (u2— 1)? a*NI/274]. (39) 


This theory is certainly not true for larger particles. It is not possible to say 
exactly at what stage the line of argument advanced in this paper cases to be 
valid, although we can take it to be applicable, if the radius is larger than 
the wave-length of light. This point will be discussed more in detail in a 
later publication. 
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I wish to express my sincere thanks to Prof. Sir C. V. Raman for 
suggesting the problem, and for the kind and encouraging guidance he gave 
during the course of the investigation. 


Summary 


The phenomenon of the transmission of light through a cloud of 
particles distributed at random in a transparent medium is theoretically 
investigated on the basis of wave-optics. The cases of both transparent 
and of opaque particles are considered, and it is found that the trans- 
mitted beam is progressively attenuated, the intensity diminishing exponen- 
tially with the increase in the thickness of the medium. The actual values of 
the attenuation coefficient are calculated in both cases. It is found that for 
an opaque particle, the diminution of intensity in the forward direction is 
actually double what would be expected from simple geometric consider- 
ations. For transparent particles, the transmitted intensity shows spectral 
variations, and this explains some of the phenomena hitherto not well 
understood, such as the colours shown by the transmitted light, and the 
complementary nature of this colour and the colour of the light scattered by 
the cloud in the forward direction. 


REFERENCES 

1. Balakrishnan, T.A.S.  .. Proc. Ind. Acad. Sci., A, 1941, 13, 188. 
2. Barus, C. .. Am. Journ. Sci., 1907, 24, 309. 

3, —————— .. Ibid., 1908, 25, 224. 

4. Bromwich, T. J. 1’A. .. Phil. Trans. Roy. Soc., 1920, A 220, 175. 
5. Debye, P. .. Ann. d. Physik., 1909, 30, 57. 

6. Mallock, A. .. Proc. Roy. Soc., 1919, A 96, 267. 

7. Mie, G. .. Ann. d. Physik., 1908, 25, 377. 

8. Rayleigh, Lord .. Phil. Mag., 1899, 47, 375. 

9, ———__— .. Proc. Roy. Soc., 1910, A 84, 25. 
10. Richardson, L. F. .. Ibid., 1919, A 96, 31. 

































TREND ANALYSER 


K. C. SREEDHARAN PILLAI 
(Travancore University) 


Received March 23, 1943 
(Communicated by Prof. K. B. Madhava, F.A.sc.) 


1. Introduction 


Tue Analysis of Time Series is of great importance in economic and other 
social studies. The eye of the economist always rests upon the Trend of 
prices and closely follows its rise and fall. The sales manager takes keen 
interest to know how the volume of sales is changing, when and why it 
fluctuates. Thus the plans of business and schemes of economic forecasting 
are all based on the study of Trends. 


Several methods appropriate to the Analysis of Time Series have been 
developed. A straightforward and easy procedure is the graphical 
method. The graph may show a general trend but it lacks in precision 
and accuracy. Hence we resort to more refined methods of analysis. 
From the practical point of view the different methods of getting the trend 
may be looked upon as methods of graduation or smoothing. These 
methods may be classified under the following heads :— 


1. Method of Moving Averages. 

2. Method of Differences. 

3. Method of Least Squares. 

4. Method of Orthogonal Polynomials. 


We discuss these methods with the help of a numerical example and 
finally give a novel method of getting the trend by a mechanical device. 


2. Methods of Graduation 


Column 2 of Table I gives the figures of Bank Clearings in New York 
City for 41 years from 1860 to 1900.1 These figures have been graduated 
by the different methods and the results are given in the remaining columns. 


In the Method of Moving Averages the average value for a consecutive 
number of years is taken as the normal or trend value for the middle of the 





1 Frederic C. Mills. Table 59, p. 259. 
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period covered in the calculation of the average. Columns 3 and 4 give the 
results of the five-year and nine-year moving averages. Fig. | shows 
the trend as determined by the Method of Moving Averages. The nine- 
year average curve is smoother than the five-year one. Either of these gives 
more or less an idea of the general trend. 








) 1 


t oe a 
8 88 90 92 94 96 98 1900 





4 i t Fy i t 1 Ain 
70 72 74 76 78 80 82% 


Fic. 1 
Observed values (Bank clearings). 
—~_—~. 5-year moving average. 
-.-.--— 9-year moving average. 

It is however desirable to get the trend based on a more mathematical 
basis.2. Attempts in this direction have given rise to several methods which 
involve the successive differences of the observed values. We shall consider 
only one of these—the Summation formula of Spencer—namely, 


¥o= [5] - [5] - [7] - (1 — 48) + ug = 125. 
The graduated values by this formula are given in column 5. 
The methods of graduation based on differences have been studied in 
detail by Whittaker and Robinson.* They offer the following remarks 


regarding these.—‘* The methods described above fulfil their purpose of 
smoothing out irregularities from observational data, in a way which is on 


2 The one fundamental difference between moving averages as a measure of trend and the 
mathematical curve is that the one implies no definite law and is a flexible measure of trend adopting 
itself to changing conditions, purporting to be nothing more than an empirical approximation 
to the drift of the series, while the other assumes a law of change underlying the variations and 
postulates a uniform and consistent trend capable of mathematical expression. 


3 Whittaker and Robinson, p. 303. 
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the whole efficient. From purely theoretical standpoint they are 1 ot alto- 
gether satisfactory since each of them contains arbitrary or empirical 
elements whose introduction does not appear logically necessary.”’ These 
remarks are true of Spencer’s formula. Another difficulty about the 
method is that we get only a small number of graduated values after all the 
labour in numerical computation. 


We may now proceed to discuss the method of Least Squares which, 
of course, rests on a logica! basis. This method may be briefly stated thus : 
“ Supposing the ungraduated values u, to be plotted as points against the 
corresponding values of x, fit a parabolic curve of some assigned degree j 
to the points (U_y, U1, Ug)----U_) determining the constants of the curve 
by the method of Least Squares aid take the ordinate of this curve x= 0 
as the graduated value of uy.’* If j be equal to 2k or 2k + 1, we can fit 
a straight line to the data taking k =0. Putting kK =1, we can fit a 
parabola of degree 2 or 3, and similarly for other values of k. To graduate 
the data at hand, we take x = 10,k = 2. The graduated values are given 
in column 6. The numerical work involved in this method is not anything 
less, if not more, than in the method of Spencer. We see also that the 
end values could not be graduated by this method. 


In the methods discussed above, the trend is not given as a mathematical 
function. It is possible to derive a mathematical equation to represent the 
trend by assuming that it is given by 


Y =a+bx+ex*?+d°+.... 
It is however advantageous to take the equation in the form 
Y =A + Bu, t+ Cus+ Dugt.... 
where 4, flg,.... are a set of orthogonal polynomials. This modification 
has been systematised for numerical work by Fisher. The several approxi- 


mating curves are obtained in successive stages, the success of fitting terms 
of higher and higher degree being tested for significance at each stage. 


The data given in Table I are graduated by this method and the 
graduated values for the Third and Fourth degree parabolas are given in 
columns 7 and 8. 





4 Whittaker and Robinson, p. 291. 
5 Fisher, R. A., p. 140. 
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I TABLE I 
Bank Clearings of New York City and Graduated Values 
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32-55 
33-79 
35-46 
38-78 | 37-61 
40-84 | 40-29 
43-21 | 43-58 
45-91 | 47-53 
48-97 | 52-21 
































3. The Equilibrium Method of Schumann 


T. E. W. Schumann has recently devised a new method of smoothing 
which he calls the equilibrium method. His method is as follows :— 
Imagine an inextensible string tied to the points u, and u, (the crude data 
being y:, Vay---- Yn) the corresponding graduated values of y, and y». 
Initially the string is not stretched tightly but is supposed to exhibit a 





*T. E, W. Schumann, Phil. Mag., 26, 179, pp. 970-83, 
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certain degree of slackness. Imagine springs to be attached at the points 
y, tO y, and to the string, each spring drawing the string in a positive or 
negative y direction with a force proportional to its own length. The 
equilibrium position of the string is the smooth curve for the data. By 
virtue of the equilibrium position Schumann equates the forces to zero and 
forms an equation which after some substitutions is brought to the form 


ap >= oa [ y, + ria + Yr) a/p + (Vrse + Yr-2) a*/p? aut ee mF 


where a lies between 0 and 0-5 and p is given by the equation p?— p + a?= 0. 
Schumann gives a few tables of coefficients to get the graduated values. 
Using the tables the data in Table I for 15 years from 1874 to 1888 has 
been graduated. The values obtained are given in column 9. (a = 6/13, 
p = 9/13.) 

Schumann admits that his formula does not apply to the ends of the 
series. However ‘ the method exhibits some rather attractive features such 
as comparative ease of calculation, and the absence of any arbitrary manipu- 
lation to include the end terms in the graduation.’ He also shows that 
close relationship exists between the method of Whittaker and the equili- 
brium method. ‘It has an added interest in that it gives a clearer insight 
into what the smoothing process actually involves. In the second place 
it points the way towards the construction of a device whereby the smooth- 
ing of time series might be done mechanically without any lengthy calcula- 
tions. It is not intended to describe such a mechanical device, but 
theoretically at least, it is possible to construct one.’ 


4. The Mechanical Analyser 

The work of Schumann has suggested the possibility of constructing 
a mechanism that may be used to smooth a time series. This has made 
practical what Schumann feared would only remain theoretical. In 
principle this device consists of a number of smooth rings that may be 
adjusted to occupy positions corresponding to the points plotted from the 
data. This adjustment is effected by mounting each ring to a rod carrying 
a spring at either end. By this arrangement each ring can be made to take 
any desired position by stretching the springs and further about this posi- 
tion, it is capable of undergoing small oscillations along the line of the rod. 
A steel wire is now passed through the rings and the curve along which the 
wire comes to equilibrium gives the trend. 

The mechanism’ has been constructed in the Department of Mathe- 
matics of the University and is able to graduate 50 values spaced at equal 





? This device has been constructed on the suggestion of U. S, Nair, 
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intervals (Fig. 2). It consists of a wooden frame 5’ x 3’ mounted vertically 


and from the top are suspended 50 springs each 10” long at equal intervals 
of half an inch. To each spring is connected a brass rod 16” long which 
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moves vertically between holes drilled in two thin metal bars fixed hori- 
zontally to the frame 8” apart. At the middle of each of these rods is 
mounted a small ring. The lower extremity of the rod is connected to 
another spring 5” long and this spring carries a rod which may be pulled 
vertically and fixed in any position by means of screws. A movable plat- 
form to hold the graph paper on which is plotted the observed points fits 
into the space between the horizontal bars. 


The working of the apparatus is simple. The observed values are 
plotted on a graph paper (scale 4 inch as unit). The graph paper is fixed 
to the platform so that the points lie along the vertical rods. The rods are 
adjusted so that the rings take the positions of the points. A thin steel 
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wire is passed through the rings. The wire is then given a slight tension 
by pulling at both ends and allowed to come to the position of equilibrium. 
The new positions of the rings are plotted on the graph paper. 


To make allowance for the elastic properties of the wire, we assume 
that the trend values are proportional to those as determined by the wire 
and calculate the constant by minimising the X* obtained. Thus, if y,; and z; 
give the observed and the corresponding wire values, 

x= 2 (y; — a + 2;)%/az; 
Minimising this expression, 
a*= 2 (y,;"/z,) + (22;). 
Applying this method to data in Table I, the values in column 10 are 


obtained. Fig. 3 shows the graph of the trend as determined by the Trend 
Analyser. 
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———— Observed values (Bank clearings). 
Values by Trend Analyser. 














We may remark that it is possible to obtain any degree of smoothness 
with the Trend Analyser by using thicker and more inelastic wires. Thus, 
we are able to pass from the straight line to more and more complex 
curves to represent the trend. Also the actual curve can be drawn by 
passing the pencil point along the equilibrium position of the wire and finally 
shifting the axis of the argument through a distance corresponding to the 
value of a. 
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3. Conclusion 


In conclusion, the different methods of smoothing have been compared 
with the help of a numerical illustration and a novel method devised to 
determine the trend. The mechanical device has several advantages, the 
most important of these being its simplicity. Also without any lengthy 
calculations, it gives readily a series of curves having varying degrees of 
smoothness. 


The author wishes to thank the University of Travancore for allowing 
him to carry on this work. He is also grateful to U. S. Nair for valuable 
advice he has given in the preparation of this paper. 
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Gnetum scandens (Malayalam—QOdal) is a lofty dioecious climber found 
commonly in the hilly tracts of the Western Ghats. The oil expressed from 
the seeds is used, in Travancore, for massage in cases of rheumatism, for 
lighting and, to a small extent, as an edible oil. 


The seeds used in this work were obtained from North Travancore. 
They are oval and yellowish with a thick shell (70%) and a kernel (30%), 
which is firm and yellowish-brown in colour. 


Experimental 


Finely powdered kernels, dried at 105° (86 gm.), were successively 
extracted in a soxhlet apparatus with various solvents (Table 1). 


All the extracts were again extracted with light petrol (b.p. 50°-60° C.) 
and then with alcohol and thus resolved into an oil and a reddish-violet, 
non-crystalline solid. 











TABLE I 
Order of Name of solvent Weight in grams | Percentage extracted 
extraction 
1 Light petrol (50°-60° C.) 12-3 14-2 
2 Ether “A 2°8 3°4 
2 Chloroform 3-7 4-3 
4 Ethyl acetate 2-7 3-2 
5 Alcohol 14-0 16-2 
35-5 41-3 














The ash of the kernel contained carbonate, phosphate, sulphate, iron and 
potassium. 


The oil had the following constants :— 


Sp. gravity 30° C. ta . "7 0-9251 
Ref. Index 30° C. be a oe 1-4694 
Acid value ee we oe ee 
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Saponification value " be -- 198-9 
Iodine value (Winkler) a .  —s 
Hehner value or es -. eax 
Acetyl value ‘in ms -- 26°98 
Unsaponfiable matter = oa ia 0-81% 


The mixed acids—These were prepared by hydrolysing about 50 gm. 
of the oil with alcoholic potash and acidifying the soap so obtained, after 
removal of the unsaponifiable matter. A semi-solid mass (mean mol.wt.= 
280-1), liquefying completely at 32°C., was obtained. 


This was resolved into solid (70%) and liquid acids (30%) by Twitchells’ 
lead-salt method. 


The liquid acids had the following constants:—Mean mol. wt =281-4, 
Ref. index (30° C.)= 1-4637, and Iodine value= 110. 


Bromination of liquid acids—A well cooled solution of bromine in 
glacial acetic acid was added to an ethereal solution of the liquid acid (5-66 
gm.) containing a little glacial acetic acid till a permanent colour remained. 
No insoluble hexabromide separated on keeping over ice. The ether was 
removed and the oily bromide was dissolved in hot light petrol (50°-60° C.) 
On cooling in ice, a brown powder (1-04 gm.) separated (m. pt. 112°- 
113° C.; after re-crystallisation from petrol it melted at 113°-114°C.). This 
is therefore linoleic acid tetrabromide. Bromine content (Stepanow’s 
method) found =52-8%; calculated = 53%. 


The filtrate was freed from light petrol, washed with dilute solution of 
sodium thiosulphate to remove any excess of bromine and then again dis- 
solved in light petrol and dried. The petrol was removed and a liquid 
bromide resulted. This might be oleic dibromide. 


Oxidation of liquid acids.—3 gms. of the liquid acids were dissolved in 
dilute alkali, the volume made up to 300 c.c. and a concentrated solution of - 
1-5 gm. of potassium permanganate solution added slowly, with continuous 
stirring, in the course of 15 minutes. The excess of permanganate was 
destroyed by passing sulphur dioxide. A white precipitate was formed. 
This was filtered, washed with a little ether to remove any free liquid acid 
and dried. The dried solid was extracted with cold ether to remove the 
tetra-hydroxy acid. The residue melted at 131°C. and had a_ molecular 
weight of 314-3. Hence it must be di-hydroxy-stearic acid. (Molecular 
wt.= 316.) The presence of oleic acid is thus confirmed. 


The solid acids.—The solid acids (13 gm.) were separated into fractions 
by fractional crystallisation from alcohol and water. 
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TABLE II 
Fraction Weight Melting point Molecular weight 
gm. "Cc. grms. 
1 5-0 68-69 283°5 
2 4°11 62 281 
3 ocd 1-49 
4 Residue 2-40 

















Fraction 1 consists of stearic acid m. pt.=68°-69°C.; Mixed melting 
point —68°-69° C. 


Fraction 2 was separated fractionally into four parts. 








TABLE IIT 
Fractions Weight Melting point Mean mol. weight 
gm. gs gms. 
A 0-811 69 282-7 
B ee “ 5s 0-362 68-69 280-3 
C ay e ee 1-697 35-36 
D Residue se ae 1-240 30 














Fractions A and B (Table III) were confirmed to be stearic acid. Fractions 
C and D (Table III) were mixed with fractions 3 and 4 (Table II), and 
again resolved by the lead salt-ether method. The solid acids so 
prepared (about 4 gms.) were fractionally separated. 








Weight Melting point Molecular weight 
gm. "GC 
x 0-95 59 264 
: - 1-21 60 257 
Z Residue 1-84 61°5 254 














‘X’ is a mixture of palmitic and stearic acids. From table* it is a 
mixture of 10% stearic acid and 90% palmitic acid. Y gave a mixed melt- 
ing point of 60°C. with palmitic acid and Z did not depress the melting 
point of palmitic acid. 

Methyl esters of the solid acids—The solid acids, prepared from a fresh 
quantity of the oil, were converted into methyl esters, in the usual manner 
with 3% methyl alcoholic hydrochloric acid. After drying in vacuum, 





* Allen’s Commercial Analysis, Vol. IV, p. 312. 
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20 gm. of the esters were distilled under reduced pressure (5 mm.) and five 
fractions were obtained. 

All these fractions were saponified separately, the corresponding acids 
were liberated and these were fractionally crystallised to isolate the indi- 
vidual acids in the usual manner (Table V). 











TABLE V 
Fractions B.P./5 mm Weight atiee Melting point 
wee mere ‘ weight of the acid 
~~. gm. *c. 
S; “ rar = 185 -190 0-5 254 60-61 
S. oe ae “4 190-195 3-0 258 61 
S; + = a 195 -200 3-5 266 58-59 
S, me ..| 200 -210 8-0 283 67-68 
Ss i i 210 -215 3-5 282 68 
Residue (by difference) os 2 a 67-68 

















The melting point of the acids after two crystallisations from dilute 
alcohol, are given in Table V. Those from S, and S, did not depress the 
melting point of pure palmitic acid and from S, and S,, of pure stearic acid. 
S, is presumably a mixture of palmitic and stearic acid. The acids from 
the residue after four crystallisations gave a product melting at 67°-68° C. 
Thus only palmitic and stearic acids were present. 


Unsaponifiable matter —This was obtained as a waxy, brownish yellow 
mass. It gave the colour reactions for phytosterols. 


Summary 


The chemical and physical constants of the fatty oil from the seeds 
of Gnetum scandens (Roxb.) have been studied. The oil consists of the 
glycerides of oleic (27%), linoleic (3%), palmitic (14%) and stearic (56%) 
acids. 

The author acknowledges his indebtedness to Dr. K. L. Moudgill for 
his interest in this work. 
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In this paper necessary conditions and sufficient conditions for that a given 
group can be ordered, will be investigated. It has been stated already in 
an earlier paper! of the author that in an ordered group the condition (a) : 
“every element = 1 is of infinite order’ is satisfied and that every group 
in which holds (5): “the commutatorgroup lies in the centre, and (a) is 
satisfied ’’ can be ordered. It will be shown that the necessary condition 
(a) is not sufficient and that the sufficient condition (b) is not a necessary 
one. Moreover some more necessary conditions will be established. 


§1. Let a, 6 (with and without indices) run over the integral numbers. 
Construct a group G of elements {a, b} with the unitelement {0, 0} by the 
following law of composition : 


{@, by} {@2, b,} = {a, b}, 
where 
a= a,+ as, 


b= (— 1)*? b, + dy. 
Then G is a group since every element {a, b} has an inverse element 


{— a, (— 1)**" 5} and the associative law holds. The product of 3 elements 
is of course independent of the bracketing: 


{@y, by} {@2, be} {a3, bs} = {y+ det as, (— 1)%** by + (— 1)*bg+ 53}. 
Moreover for every positive integer m there is 
{a, b}” = {ma, b [1+ (— 1)?+....+ (— I)" }} 
and the right hand side is equal to {0, 0} if and only ifa=b=0. Hence 
condition (a) is satisfied in G. However 
{1, 0}-? £0, 1} {1, 0} = (0, — 1} = {0, 1}. 


Suppose G is ordered and {0, 1} a positive (negative) element, then its trans- 
formed elements are also positive (negative) whereas the inverse element 





1“ Ordered groups,” Proc. Ind. Acad. Sci., 1942, 16, No. 4, Sec. A. 
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is negative (positive). Thus G cannot be ordered and therefore (a) is not a 
sufficient condition. 

§2. Let R be a subring of an ordered field with the property that 
when a belongs to R, then (a?+ a): 2 belongs to it. E.g., R may be an 
ordered field, say a field of real numbers, or it may be the ring of the integral 
numbers. Let a, b, c, d, e run over the elements of R and construct a group 
H of elements 

[a, b, c, d, e] (1) 
with the unitelement [0, 0, 0,0,0] by the help of the following law of 
composition : 

[a;, by, C1, dy, €;] [a2 be, Co, dz, 2] = [a, b, c, d, e], where 
a=4+4, b=b,+b, c=¢,+ C2— agh, 
d — d, + d2+ byd2(ad2+ 1): 2- CQ, e= ar €3+ byaob.+ ab, (6b, +1): 2-— C15, 


[a, b, c, d, ey =[— a, — b, — c— ab, — d—ba(a— 1): 2— ac, —e— ab(b—1):2— bc] 


That the associative law holds, can be checked by an elementary calcu- 


lation. Hence H is a group. Put 


[1, 0, 0, 0, 0] = a, [0, 1, 0,0, 0]= 8. Then the commutator 

(a, 8) = aBa-1p-4 = [0, 0, 1,0, 0] and 

(<, (a, B)) = [0, 0, 0, 1, 0), ((8, (a, B)) cs [0, 0, 0, 0, 1} 3 
thus (a, 8) does not lie in the centre since it is not commutative with a and 8. 
Therefore H does not satisfy the condition (b). However H can be ordered 
by the following rule : An element (1) is considered to be positive or negative 
according as the first of those of its five co-ordinates which are different 
from zero, is positive or negative. The unitelement is considered to be 
both positive and negative. 


Now (2) shows that the product of two positive elements is positive 
and (3) shows that inverse elements have opposite signs. If one transforms 
(1) by any element of H, one gets an element of the form 


[a, b, c+ fi, d+ fo, e+ fal, 
where (as can be shown by an elementary calculation) a= b=0 implies 
1=0 and a=b=c=0 implies f,=f;=0. Hence transformed elements 
have the same sign ; thus all the conditions for an ordered group are 
satisfied. Hence the sufficient condition (b) is not necessary. 
If in particular the ring R is the ring of the integral numbers, H is 
generated by a and 8. It may be mentioned without proof that in this case 
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H is isomorphic to F/C;(F) where F is a free group with 2 generators and 
C; denotes the third commutator group. 


§ 3. It has been shown already in § 1 that a group in which an element 
+1 is a transformed of its inverse element cannot be ordered. From this 
remark it follows immediately that if a manifold contains amphidromic 
curves which are not nullhomotopic, its fundamental group cannot be 
ordered. Of course the groupelement which corresponds to the amphi- 
dromic curve taken with a certain orientation is the transformed of the 
groupelement corresponding to the same curve taken with the opposite 
orientation, i.e., to its inverse element; as the curve is not nullhomotopic 
the groupelement is different from 7. Hence the fundamental group of a 
one-sided surface cannot be ordered. 


Moreover in an ordered group an equation x”=c cannot have more 
than one solution for any positive m, since if a and 5b are two different 
solutions and a> b, then a” > b”. For c =7 this condition is identical 
(3 with the condition (a). 


Another generalisation of condition (a) is the following theorem: 


Let a,, ..++, Ag be n elements of an ordered group A and not all equal 
to 1, let B be the subgroup generated by them and C(B) its commutator group ; 
then B/C(B) contains elements of infinite order. 


Proof.—Suppose that the highest of the absolute values |a,|, ...., |as| 
is equal toa. Let A, be the group of the elements with absolute values at- 
most comparable with a, and A, the subgroup of the elements with absolute 
values infinitely small to a. Then A,> B anda liesin B. In A, there 
exists a function? ,/g which is zero for all the elements of Ay, but differ- 
ent from zero for the other elements of A,. Now ,/g (bcb%c)= 
aigb + lg c — lg b — ,lgc =0. Hence the commutatorgroup of A, lies in 
A, and in particular C(B) © Apo. As lg a*=n, none of the elements 
e™, a**,...... is situated in C(B). Hence the element a C(B) of B/C (B) 
is of infinite order. 
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7. Introduction 


AMONGST the best-known phenomena of meteorological optics are the 
corone consisting of a central disc and one or more coloured concentric 
rings seen surrounding the sun or the moon, when viewed through thin 
clouds. The observed angular extension of the corone is variable, being 
generally between a degree and ten degrees. Closely related to these are 
the iridescent clouds exhibiting vivid colours, which are sometimes observed 
at distances of about 5° to 30° from the sun. The correlation between 
these two phenomena was recognized by Simpson (1912), who showed that 
iridescent clouds are nothing but parts of corone. 


A mathematical theory of the corona was put forward by Airy, who 
assumed that the water droplets in the cloud act like opaque disks, and thus 
give rise to a diffraction pattern of bright and dark rings. Obviously, 
however, One is not justified in considering the water droplets as opaque, 
since actually they are transparent. Indeed, many of the experimental facts, 
as for example those observed by Barus (1907, 1908), Mecke (1920 a), and 
Mitra (1928) are not explained by the “‘ opaque disk” theory. It is evident 
that the light transmitted through the droplets must be taken into account. 


Mecke (1920 5) has tried to give a more complete theory, taking into 
account the rays that pass through the drops also. His method consists 
in dividing the effect of the droplet into three parts, due to diffraction, trans- 
mission and reflection respectively. The first part is taken to be the same as 
that due to an opaque disk of the same size. The transmission effect is 
calculated by tracing the refracted rays, and finding their divergence and 
phase change. A similar calculation is also made for the rays reflected 
from the surface of the droplet. All these parts are taken to be coherent, 
the amplitudes are added, and the resultant intensity is calculated. 


It may be pointed out, however, that Mecke’s application of the ideas 


of geometrical optics cannot be justified when the particles under consider- 
ation have a radius which is only a moderate multiple of the wave-length. 
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For, the transmitted rays themselves form a beam whose aperture is of the 
order of the diameter of the droplet, and hence comparable with the 
wave-length of light. Such a beam would give diffraction effects of its 
own, so that at the focal plane, where the geometric rays are supposed 
to come to a focus, the diffraction pattern would have a size about the 
same as that of the droplet itself. The ray-paths indicated by geometrical 
optics cease, therefore, to have any significance when one is dealing with 
droplets whose radius is only a small multiple of the wave-length. It is 
obvious, in these circumstances, that one must treat the problem on the 
basis of wave-optics alone, and not on the lines suggested by Mecke. 


A theory of corone based exclusively on the ideas of the wave optics 
has been outlined by Balakrishnan (1941) taking into consideration the 
phase changes occurring in the passage of the light through the drops. 
The present author has used similar ideas in discussing the transmission of 
light through a cloud of transparent droplets, in a paper appearing earlier 
in these Proceedings. Both in that paper and in Balakrishnan’s work, the 
aspect considered was the evaluation of the intensity in the forward direc- 
tion, viz., along the incident rays. The present paper concerns itself with 
the evaluation of the intensity in other directions, which is necessary to 
find the position and intensity of the rings observed in the corona. 


2. Derivation of the Formule 


It has been shown by the author (loc. cit.) that the amplitude of the 
wave diffracted by a droplet of radius a in a direction making an angle 
@ with the incident direction is 


X = X, ~7 Xz, (1) 
n[2 
where X,= +) Jo (n sin 4) sin (X — € cos 6) sin @ cos 8 dé, (2) 
0 
Xp= A 2. sin y, (3) 
& = 4n(p— 1) a/A, and y= 2m (sin 9)/A, (4) 


the incident wave being represented by sinX. The whole problem thus 
reduces to the evaluation of the integral in the expréssion for X,. The inte- 
gration is not simple for finite values of 7. The author tried various 
methods of putting it in a form amenable for numerical computation. The 
best method was found to be one of partial integration, by which it could be 
put in the form of a series, of which the later terms are small compared with 
the first. The method consists in partially integrating sin (X — € cos 8), 
A4d 
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keeping J, (y sin #) as a constant. We may then write the integral as 
n[2 
= af SoC sin & sin @ — & 605 0) (6 c0s 4) d(é cos 0, which i 
5 


ria 
0 


= - [40 (n sin ) {sin (X — € cos #)— é cos 8 cos (X — é cos | 
— | involving J, 
- n [Jo (n) sin X — sin (X — €)+ & cos (X — | _ f involving J,. 


Normally, the value of é is fairly large, and we can neglect all terms 
invloving higher powers of 1/€ as compared with those involving 1/é. Then, 


the integral appearing in X, is = : cos (X — €), so that 


X= K [cosX cos (€)/€+ sin X {sin (£)/E — J, (»)/n}] (5) 
where K = 27a?/A. Hence, the intensity is 
I= K? [{cos (£)/€}? + {sin (€)/é — J, (n)/n}?] 
L , J2(m) _ 25, (y) , sin J 
oo Eth 1 ied «=F i. 
“— lz oe, ” g 6) 
This expression (6) is not very suitable when é is small, and also when 
7 becomes comparable with €. In such a case, a second approximation can 
be worked out by expanding the integral for X, and summing up the most 
important terms (vide Appendix I) which gives 


1=k+[ l J,’ () _ 97h () Z sin (€+ mits). (7) 


at 7 7 g 
A third approximation can also be obtained, which leads to the expression 
for I as 








[= K? [S*+ {C— J, (m)/n}*], 
where S= : cos (€+ 7/2) 


it ., . n* P - .. . 
+ ga| sin (6+ 01/28) —% cos (+ 1°/26)— gh sin (€ + 19/24)] 
and C= } sin (§ + 9°/2é) 
! 1 2 * 1 2 7f 2 ° 
+ gal cos (6-+ 09/24) + sin (€-+ 9°/28)— ghs cos (E+ 07/28) 


3. Discussion of the Expression for the Intensity 


For large values of €, the first approximation (6) is sufficiently accurate 
to evaluate the positions of the maxima. In order to determine how far it 
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is different from the more accurate expression (7), the author actually 
computed the intensity distribution given by both of them for a few cases. 
It was found that the higher degree of approximation did not appreciably 
alter the positions of the maxima and the minima, but that it gave values 
for the actual intensities at these positions different from those given by 
(6). Therefore the simpler expression (6) could be used for a discussion of 
the positions of the bright and dark rings in the corona, and, in the follow- 
ing discussion, this alone has been used. 


We shall first derive a few general properties of the function on the right 
hand side of (6). For drops of a particular size, € is a constant, and the 
quantity which varies with the angle of diffraction ¢ is 7. Hence, the 
varying part of I is given by 

I’ = K* [J,* ()/n? — 25, (m)/n * sin (€)/€] (8) 
Putting sin (€)/€ as k, this can be written as 
I’ = K? [J,* (m)/n? — 2kJ, (m)/n] (9) 
Differentiating with respect to 7, 
dl’ 
= 2K? [J, (9)/n— k] + Je (m)/- (10) 


Hence, the maxima and minima of I’ occur for the values of » given by 
Ji (in=k (a) 
Js (y)/n=0 (0) 
Call the roots of the former a,, d2,.... and those of the later b,, be,.... 
Differentiating (16) once again, we get 
[470 /dn*)n=am = + 2KIF27(y)/? (11) 
which is always positive. Hence, at » =a,,, I is always a minimum. 
Again, [4*I’/dn*]n2o = — K*(1 — 2k) [d/dn {Jx(m)/n}n20 
= — K*%1 — 2k)(a +ve quantity) (12) 
This is positive or negative according as k is > or < 1/2. Hence, at n = 0, 
lis a maximum or minimum according as k 2 1/2. 


The other roots of (6), namely 5,, by,.... may correspond to a maximum 
or minimum. It is also easily seen that they correspond to the value of 7 
for which J,(y)/n has a turning value. 


4. The Rings of the Corona 


We are now in a position to discuss the disposition of the rings in the 
corona, which correspond to the maxima and minima of intensity. In this 
discussion, we shall limit ourselves to the case when monochromatic light is 
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used. The nature of the ring system with white light can be readily 
deduced from this. 

The diffraction in the forward direction and the colours exhibited by 
the central portion of the corona have already been discussed in the previous 
paper. It is now only necessary to consider the outer rings of the corona. 
First of all, it can be shown that the thickness or density of the fog has no 
influence on the position of the rings in the corona. The intensity of the 
light diffracted by a single droplet at an angle ¢ is given by (6). But the 
light reaching it, and the light diffracted have both to travel through the 
thickness of the fog. This total path will not differ appreciably from the 
total thickness / of the fog, if the angle ¢ is not very large. Hence, from 
the equation (32) of the paper already cited, the intensity of the light 
diffracted per unit volume of the fog in the direction ¢ is 


i = NK* [z: A oe - re ' aad exp {- 4na°NI(, _ nt \ 
(13) 


where N is the number of droplets per unit volume. If the size of the droplet 
is fixed, € is a constant, and the position of the ring system is determined 
only by the expression within the square brackets. Neither N, nor / enters 
this expression, so that the ring system must be uninfluenced by the thickness 
of the fog. This has already been verified by Barus (1912) and also by 
Mecke (1920 a). 

Taking, therefore, the expression within the square brackets, we have 
already shown that the position of its maxima and minima are given by the 





8, 8s 








Fic. 1. Graph of J, (») / » against 


roots of the equations (a) and (5). In the following discussion we shall limit 
ourselves to the first two or three rings, since, in general, the outer ones 
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would be faint. Also, from Appendix III, it will be seen that the series 
for C and S are not convergent if » > é, so that this discussion cannot be 
extended to values of »>&. This corresponds to sin ¢ > 0-667, or 
¢ > 41°49’. Since the corone usually observed do not extend beyond 
about 20°, this does not present any serious difficulty. 

A careful analysis of the case shows that the existence and the number 
of the roots of (a) and (b) depend very much on the value of k = sin (€)/é, 
and that there are four distinct cases. Denote by B,, By,.... the maximum 
and minimum values of J,(y)/n which correspond to 7 = b,, by, bg, . 
The course of the curve J,()/n is shown in Fig. 1, and the values of the B’s 
and the b’s are tabulated in Table I below. 


TABLE [| 





b= 5:14 b,= 8-42 bs= 11°62 
B,= —0-0662 B,=+0-0322 =— 0-200 


by= 14-78 b;= 17:96 
B,=+0-0140 B,=— 0-0105 

















The four cases are as follow: 

Case I: k is +ve, and 1/2>k>B,. 

In this case, equation (a) has only one root, a,, between zero and 5b,, 
which corresponds to a minimum value of I. Consequently the distribution 
of maxima and minima are as shown in Table II (a). 

Case Il: kis + ve, and k < By. 

Equation (a) will have 3 roots at 7 =a, ad, and a3, where 
a,< b,, and a, and a, are on either side of b,. The arrangement of the 
maxima and minima is tabulated in Table II (0). 

Case III: k is —ve, and | k|< |B,| 

In this case, equation (a) has two roots, a, and az, on either side of b,, 
which give rise to the arrangement listed in Table II (c). 











TABLE II 
Max 0 b, bs 
(a) Case I 
Min. a bg bs 
Max. 0 b i by! b 
(b) Case 1 wee te oe 3 
Min. a, \ a a, | be 
Max 0 b, | bs by 
(QCacII .. witli —— 








(d) Case IV 
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Case IV: k is — ve, and |k| > |B,| 
Here, (a) has no roots at all, and the maxima and minima occur at 
n = b,, bg,.... as shown in II (d). 


The dotted lines are intended to show that the maxima within them 
from a sort of gradual transition from the case before to the next case. This 
point will become clear if the intersection points of the graph y =k and 
y =J,(y)/n are studied. In case I, there is no intersection point beyond 5,; 
but as k decreases, case II becomes operative. At first, a, and ag are close 
together, so that the maximum at 5, will not be prominent. The reverse 
happens in case III, the first maximum becoming less and less pronounced, 
until it becomes a minimum in case IV. 


It is thus clear that the nature of the ring system depends very much 
on the value of k, so that even a small fluctuation in the radius of the droplets 
produces a large change in the disposition of the rings. It is therefore 
interesting to see how the rings behave as the particle size steadily increases. 
Suppose that the wavelength A of the light is a constant. Let us begin with 
the value a/A = 2-25 for which € = 3”. This is a case midway between 
cases II and III. There will be maxima at y = 4,, by, etc., and the rings 
will appear to be close together. As & increases, k becomes negative, case 
III becomes operative, and the first ring becomes fainter and fainter until 
it disappears for |k| > |B,|.. The whole ring system will now appear to 
have expanded. This continues for some time, after which k numerically 
decreases, and the changes take place in the reverse direction, the rings 
coming close together at § =4. Thereafter k is +ve, and the second 
ring becomes more and more indistinct until it disappears when k is > B,. 
This continues for some time, after which the second ring reappears, the initial 
stage being reached once more at £ = 5z. This finishes one cycle of changes. 
Another cycle, exactly similar to the above takes place in the range of values 
5m to 7x of &, except that, in this cycle, case IV cannot be realised at all, 
since |sin (€)/€| is never > |B,|. This continues upto é = 9x, after which 
case I cannot occur. For droplets of larger radius, i.e., a/A > 7, the expand- 
ing and contracting of the rings will become less and less noticeable, until, 


for very large drops, the ring system will correspond to the case of an opaque 
disk. 


In this limiting case of very large drops, 1/€* and sin (€)/€ become 
intrinsically small, so that (10) reduces to 

[ = K* J,*(m)/n* (14) 

which is the same as the expression for an opaque disk. An interesting 

point arises in this connection. As shown by equation (34) of the previous 
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paper by the author, the decrease in the forward intensity due to a large 
droplet is the same as that for an opaque disk of the same size. Thus, for 
large droplets, the total energy in the corona must be equal to 2za*, or 
double the amount in the corona due to an opaque disk. However, the 
expression (14) shows that the corona is identical with that given by an 
opaque disk of the same size. The discrepancy can be explained by the 
reason that our theory is not valid for large angles of diffraction. The 
theory only predicts that for small angles of diffraction, the corona is identical 
with that produced by an opaque disk. There should be an equal amount of 
energy diffracted at larger angles, of which no account is taken by the theory. 


5. Numerical Computation 
In order to bring out the facts discussed above in a clear manner, the 
intensity distribution in the corona has been calculated for a few cases, using 
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the more accurate expression (7). The wave-length of light was assumed 
to be 5000A.U., and the computation was done for five values of the 
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radius of the drop. In Fig. 2, the results have been plotted with the angle 
of diffraction as abscisse, and the intensity as ordinate. The curves clearly 
show how the nature of the ring system is altered by changes in the size of 
drops. They also exhibit the contraction and expansion of the rings as 
the size of the droplets is increased. 


6. Comparison with Experiment 


The theory developed in this paper explains most of the observations 
made by previous workers in the field. Even very early measurements, such 
as those of Kamtz (Pernter, 1922) showed that the opaque disk theory 
is only asymptotic in its application, being more and more correct, the 
larger the size of the drop, as is to be expected from the theory. Mecke - 
(1920 a) and Mitra (1928) have observed that the colour sequence in the 
outer rings, with a white source, had no special order, and that it was very 
different even with droplets of slight difference in radii. This is easily under- 
stood, since our theory shows that, with monochromatic light, the distri- 
bution of the rings is highly susceptible to even small changes in a/A. The 
case with monochromatic light has also been verified by Mitra, who found 
marked changes in the position of the rings as the size of the drop changes, 


Mecke has given some values of the relative radii of the rings in 
different cases. If ¢,, d., ¢3 are the angular radii of the first, second and 
third dark rings in monochromatic light, then sin ¢,: sin ¢,: sin 43 will give 
the relative values of » for these. They have been tabulated in Table IV 
of Mecke’s paper (1920a). An examination of it will show that the ratio, 
sin ¢,/sin ¢,, call it p, varies from 1-8 to nearly 2-4, but that mostly, it is 
in the neig!:\bourhood of 2-2. The value of sin ¢,/sin 4,, call it po, lies 

etween 2-9 and 3-4. Mecke has also found that p, varies periodically 

with the radius of the droplet, and that it becomes less than 1-83, the 
normal value given by the opaque disk theory, only for a small range 
in each cycle. All these are easily explained by our theory. 

Taking one cycle, we start with €= 3a, for which k= 0, and p, = 1-83, 
and p,= 2-8, which are the normal values. As € increases, p, will slightly 
diminish and p, will increase, as may be seen from Fig. 1. This will occur 
only for a short while, for the maximum at 5, will soon be too faint, and the 
ratios will become 2-3 and 3-6 respectively. This holds for case IV also. 
If we take case II, the ratios will be in the neighbourhood of 2-2 and 3-0, 
and on passing to case I, they become larger than this. This explains the 
range of values of p, and p, as given by Mecke’s experiments as also the 
fact that p, is seldom less than 1-83, since it occurs only for a short range 
after the transition from case II to case III. 
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It is with great pleasure that I take this opportunity of expressing my 
gratitude to Prof. Sir C. V. Raman for suggesting the problem, and for 
the keen and continued interest that he took in the investigation. 


Summary 


The commonly accepted theory of corone, based on the idea that water 
droplets act as opaque disks, is not only theoretically unsound, but is not 
also in accord with experimental facts. The theory due to Mecke is also not 
satisfactory as it based partly on geometrical optics and partly on the theory 
of diffraction. In this paper, a new theory of the phenomenon is deve- 
loped, using only the principles of wave-optics, and taking into considera- 
tion the portions of the wave-front transmitted through the droplets. The 
integrals so obtained are integrated by a suitable method, and expressions 
are obtained for the intensity distribution in the corona. A discussion of 
these expressions shows that the theory satisfactorily explains most of the 
phenomena exhibited by corone, such as the extreme susceptibility of the 
ring system to even small changes in the radius of the drops, and the 
oscillation of the ring system as the radius steadily increases. 
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APPENDIX I 
Evaluation of the Integral in the expression for X, 


The integral to be evaluated is 
[2 
A= J (y sin 4) sin (X-€ cos @) sin @ cos 0 dé 


The integration is to be done in such a way that the resulting function is 
readily adapted for numerical computation when é is large, i.e., when terms 
in 1/&*, 1/&%, etc., can be neglected compared to those in 1/€. 


For the purpose of integration, it is convenient to split the integral 
into two, and write 
n|2 


A= sin X f Jo (7 sin 9) cos (€ cos 4) sin 8 cos 8 dé 


0 
n[2 


— cos X f J (7 sin 4) sin (€ cos 6) sin 8 cos 6d6 
0 
=CsinX —Scosx (say) 
Both C and S can be expanded in the form of a series. 
Put € cos 9= x, and 7 sin 0= y. 


dy _ 46 COs ae = 
Then, > +. -23 2). 
Also, the limits of Stee are given by 


6= 0, ie. x= & y=0 
O= 2/2, i.e.. x= 0, y=». 


g 
Hence, C= ef Jo (v) cos x x dx 


Partially integrating with respect to x, 
x= 


g 
=F ra | (Cos x+xsin x) J, (| i f (eos x+ x sin x) J, (y) dy 
x=9 0 
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Substituting for dy in terms of dx, 


é 
= C= | (cos €+ ésin &)—J, @)|- af (cos x+ x sin x) + x: J, (y)/y.dx 


Again, partially integrating with respect to x, the second term becomes 


2 
—— a [o cos x+ 3x sin x— x* cos x): J, oy]: 
0 


¢ 
on BS (3 cos x+ 3x sin x— x? cos x) - 2 en 


1 
—— this is 


- [46 cos + 3é sin — €? cos j— 220) 


pH 

é 

. é 
+e (3 cos x+ 3x sin x—x# cos x) + x Jy (y)/p* + dx 


Continuing this process, we can finally put C in the form 


C= (cos € + ésin £) 
— Je (Bcos +3 Esin  — Ecos £) 


sate (1S cos € + 15 sin € — 6€? cos € — &% sin &) 





+(—1) Fp ie (ap pcos € + dp p1 €sin &é 
+ Ay p-2 §* Cos E— ----) 
_ 3a? Wn) Vn Tn) _ 
ge 9 e 
+(— Pay (2) a(7) es 


where dy 4 = = 7G ; eo to ; (Vide Appendix II) 
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In an exactly similar manner, S can also be put in the form of a series 
as below: 
— 1)P-1 2-2 


D et (ote (dp, p sin € — ay, pa Cos € — dy p-e F* sin € 


r 


e=1 


+ dp, ps &®cosE +... .). 


The above two series give the complete expansion of the integrals C and S. 
Ordinarily, for the range of particle sizes in which we are interested, € is 
large, and only terms in 1/€ need be taken. Then, 


A =sinX - sin €/€ + cosX - cos €/é, so that 
T = K*{1/é*+ J,%q)/n? — 25,()/n - sin €/€]. 


This first approximation is not sufficiently accurate if 7 is comparable 
with €. Then, the last terms inside the brackets in the right hand side of C 
must be taken into account. C then becomes 


="ll-am ht 9 ge] 


Similarly, to a second approximation, S - 008 (é +32) so that 


pa wef, + BD Ho. sin(é 941207, 
- n° 0 g 

If a further approximation were required, it can be obtained by taking 

all the penultimate terms inside the brackets in the right hand side of C and S. 

They can also be summed, and lead to the result that, to a third approxi- 
mation, 

l 

uy 

[cos (g + 99/28) + % sin (€ + 9/26) — gh cos (& + 1/24) | and 


* sin (€ + 7/2€) 


+ & 


S = 7c0s (€ + 7*/2£) 
+ i [ sin (€ + 99/26) — "cos (€ + 99/26) — gee sin (€ + n*/28) | 
Any desired degree of approximation can thus be obtained by taking 


further terms in the expansion, but one more than the second degree will 
not be needed practically. 
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APPENDIX II 
To prove that ag ~=(p + 4)!/2%(q)!(p —q@)! 


The quantities a, , depend on the coefficients of the terms in the expan- 


g é 
sion of the definite integrals, J x? cos x dx, and / x’ sinx dx. If we call 
0 0 


these respectively C, and S,, then 
Cy = Do é sin € + p, 1 cos € — p, &-* sin é—.... 
Sp = — Po & cos é + p, "1 sin € + p, &* cos € —.... 
where two terms are +ve, the next two are —ve and so on. The series 


on the right hand side has (p + 1) terms each, the last coefficient being py. 
The values of the coefficients are easily seen to be 


Po =(P)'/(p — 9)! 
It is to be noted that py = pp. 


From the method of partial integration used in Appendix I to determine 
the series for C, it can be verified that the laws of formation of the numbers 
dy, are given by the following: 


4p,0 = Po %-1, 0 
a.1 = (p — 1)o Ap1.1 + P1 4p-1,0 
Qp,2 =(P — 2)o4ps,2 +(P — 1). 4p1,1 + Pe dpa, o 


Ap, p-r= 19 Apa, g-1 + 21 Opa, po +--+» + Ppa Gps, o 

and dy, » = dp, p+: 

From these substituting the values of p,, we get, 

(1) ag 9 = 4-1, o, SO that a, y= 1, since a, 9= |. 

(2) apg, + = Qp-1, itp ag-1, oS, = Qp-1, it ?p. Also, ay, ees a SO that 
aps, t= En. 

a 
(3) a, 2 = Ap4.2+(Pp — 1) ao41 +p(p — Wags 
o~1 

= Ap4,2+(p — 1) < n+ p(p—1) 
=Agi2et+(p—1) Xn 


Hence, dy 2 = ds, 2 +3 (m _ 1) En. 
3 1 
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But, de2 =4e,1 ==n =(2 — dn, so that 
1 
Ap, 2 =E(m — 1) En. 
1 1 


(4) ay, 3= ps, s+ (D— 2) Apa, 2+ (p— 1) (p— 2) p41 
+ p(p— 1) (p— 2) 45.4, 6 


-—-1 m —1 
=ipest(e—2 2 (e—?t z n+(p—1)(p—2) =n 
+ p(p— 1) (p—2) 
= Op s+ (p— 22 (m— 1) En. 
z 1 


Hence, dy, 3 = 43, 3+ S(r- 2) z (m— 1) z n. 
4 1 1 


But, d3.3  =432 =2 (m— zn = (3-— 2) 3 (m— En, so that 
ap, =F (r— 2) Z (m— Zn. 

Proceeding in this way, it is easily seen that 
dp, g= E(t q+ 2 (s—q+ 2)-+(r— 2) E(m—1) En, 


The expression on the right hand side can be summed up by using the 
well-known summation: 





z (a+ x) (a+ x4+ 1)°--(@4+x4+r—-—D= u 


(+1) (a+ n) (a+ n+ 1)-:: 


(a+n+r). 
We then get, Dn=m (m+ 1)/2=(m+ 1)! /2(m—1)! 


E(m—1) Zn =F (m— 1) m (m+ 1)/2= (r+ 21/22 (2)! (r—2)}, ete. 


Thus, 4 9=1; a, 1=(p+ 1)! /2)-(p—1)!; 

Ap, 2=(p— 2)! /2? (2)!(p— 2)! 5 ay, s=(pt 3)! /27 (3) !(p— 3)! 
and so on, so that 

Ap, = (p+ 9)! /2% (q)!(p—4)! 
It is to be noted that this relation automatically satisfies the equation 
24, p= 4p, p-1- 
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APPENDIX III 
A Note on the Convergency of the Series for C and S 
It can be shown that both the integrals C and S, expressed in the form 


of the series we have chosen, are convergent if 7 <€. Taking first C, we can 
express it as follows: 


C= sin € {4, 9/§ — (de, 1 9°/2+ as, 9 a*/2? (2) H/F 
+ (As, 2 94/2? (2)j! + aa, 1 79/23 (3)! + ag, » 08/24 (4) ! )/é 
ste in de deh i a eat a tan Se a doe } 

+ cos € {(@,, 1+ Ge, 9 n?/2)/€*— (ae, 27/2 + az, 1 04/2? (2)! 
+ 4 9 n°/2° (3) ! )/é 


J ay 
- A {ho)- 2,273 JO + 53 f 2) _ anes +} 


= Asin + B cos £— J/£* (say) 





The series S can also be expressed similarly as 
S=Acos é—B sin é. 


It will now be shown that A, B and J are convergent if 7< €. Taking 
A, if ug be its p™ term, then 


=. [feeee ??/2? (p)!+ dpso,9—1 02+ 2/29+1 (p+ I)!+-- ] 
* Lay p—1 01°F 2/2?-* (p— I) E+ p41, p-2 778/27 (p)!+ °°: 


which can be put in the form 


Pp (p+4), p(p+4) 1 , p(p+4) (p—1), xf 
wf p? “per (3)! 1" p@F VD) (PED) ” G)it a 


Up+4 
Up 





























“ova | 7 
uz a | 14?! (p— 1) (p—2) 
at ao ee | 
Hence, Lt |“#+*|—7 If n<é, this is< 1, and the series is con- 
p->co Up rs 








vergent. A similar test shows that for B also the same cnodition holds 
good. 


Taking J, if v, be its pth term, then 


Up+1 
Up 


— 7 , 4p+1, p+1 Jp (n) , 








e Qp9 Sp-1 () 
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Substituting for the a’s, and also using the relation that, when p is large, 
Sp41 ()/J (n) = 9/2 (p+ 1), we get 


psi} 2p+1 | 7? 
Up 2p e 


Vp 4. 2 : ‘ 
Hence, Lt |“2+1/— 7 so that J is convergent if 7 < é. 
p->co Up - 


Thus, the expansion of the integral for X, in the form of the serics 
C and S is convergent, and hence valid, if 7 < é. 
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THE Bessel method of summation, or the (J,, 4) process of summability was 
first introduced by Dr. S. Minakshisundaram® as a means of summing up 
series of eigenfunctions in boundary-value problems. For the purpose 
which he then had in view, it was sufficient to prove the fundamental 
theorems of this process; and the possibility of interrelating it with other 
known processes remained open. At his instance, I undertook to study 
that problem, the solution of which is the object of this paper.* 

It is proved here that the J, process is roughly of the same scope as 
Cesaro’s process. This fact, however, neither vitiates its interest as a 
generalisation of Riemann’s and Lebesgue’s methods, nor diminishes its 
convenience for the treatment of special series. 

I owe a great deal to Dr. Minakshisundaram in the way of help, 
encouragement and criticism. 


§1. Let J, (2) be the Bessel-function of order p, so that 


t# 1? 
uO = TGF) {! ~ 2 Qt * Ta OptD pth 


i4 





Let 
a, ()=2* P(e +1) a. 


Let (A,) be a strictly increasing sequence of positive numbers tending to 
infinity. 
Definition—A series Xa, is said to be summable (J,*, ») to the sum S, 
if the series 
co 
& ag [oy On OFF 


converges for real t > 0, and integral K > 0, and if 
Lim ¢,x (4) = 5, exists, 
t>0 





* This paper formed part of a thesis approved for the M.Sc. degree of the University of 
Madras, 
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AS 
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where, 
bux (= 2 dy [Oy Oe OF 
and 
gua (= dy (0. 
When K= 1, p=}, A, = 72, the definition reads as follows :— 
A series & ay is said to be summable (J;, n) to the sum S, if the series 


e sin nt 
Ean (AT) 
an nt 


n=O 


converges for real t > 0, and if 


; sin nt 
Lim 2 a, {-—— }=5, 
t>o nt 


which is Lebesgue’s summability. 
When K = 2, p= 4, A, = 2, the definition reads thus: 
A series Z a, is summable (J,*, n) to the sum S, if the series 


. sin nt\2 
= @ ( 
™ nt 


n=0 





converges for real t > 0, and if 
, sin nt\? 
Lim Za, {——— } = s, 
t>e nt 


which is Riemann’s summability. 


When p= 4, A,=n, we obtain the generalized Riemann Summability 
considered by B. Kuttner® and M. Vignaux.® 


The case when K = 1, A, =n is important for our purpose, inasmuch 
as it is closely connected with Cesaro’s method of summation, I call this 
summability J,. 


§2. Minakshisundaram has proved, among others, the following funda- 
mental theorems: 


(2:1). If a, is summable J, then a, =o (n#+4). 
(2:2). If 2% a, is summable J,, and the convergence of the series 


oo 
2 dy, a, (nt) is such as to permit integration term by term, 


then it is summable J,, for v >, to the same sum. 


I prove the following simple theorem: 











A) 


n, 
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(2:3). A necessary and sufficient condition that = a, which is summable 
J,, Should be summable J,,_, is that the series 


es n* a, Ge+4 (nt) 
converges for t > 0, and the sum-function tends to zero with t. 
For, 


>, (t)— 27u ( (u yz ae a, +4 (nt) 


= $, (t)— 24-1 Pw) PZ a, ae n® 


= by = BPE ay BLO 


MDW) Gat [7 I, (nt) — Jy 1 (0 ] 
=2a, eo). 20-1 Py) 
= $y -1 (0). 


§3. The relation between the J, process and the process of Cesaro 
is established by the following: 


THEOREM J. Jf a series X a, is summable (C, r), then it is summable 
J,, for » > r+4, to the same sum. 


Proof. Without loss of generality, the sum of the series may be 
assumed to be zero. 


Let S (x)= = a,, forn<gx<an+l; 
v=0 


e 


and M0= 1 f 1) "=1§ (1) dt; 


S (x)= S (x); S= z a,= S (m). 


The Cesaro-summability of order r of the series X a, is equivalent to the 


statement that 
S™ (x) = 0 (x’). 


By Abel’s lemma on partial summation, 


E ay a, (nt)=Z S, [ay (vt)— a, (VHT 1))+ Sp a, (mt). (3-1) 
Now, . 5 
4 
a, (vt)— a, (v+1)=— ie (xt)- dx, 


ba 
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so that 
V+1 


S, [e, () — a, (v+ 1 0)=— fse (x): < a, (xt)> dx 


Hence from (3-1) 
v+i1 


m m—1 c d ; 
* a, 4, (nt) =— ca f $0) (x)- Te tH (xt) dx+S,, a, (mt) (3-2) 
p s 
: - ; os Oe ; . 
= Sy, 4, (mt) — f S$ (x)- dx tH (xt) dx. (3-3) 
7 


Integrating the second expression in (3-3) by parts A+ 1 times, where 
h= [r], we obtain, 


m 


n d 
2X Ay a, (nt) = S,, a, (mt)— §M ; 


Ly & a sic 
en Sp on (mt) + S¢ ) nit Oy (mt) + 


0 
hed 
. d 
+ (+t SE + (mm) 5 My (mt) 


+ (~yr* f S*+1(x) D,f+3 a, (xt) dx. (3-4) 
° 
By hypothesis, we have, 
S..= 0 (n’) 
S‘®) (m)— 0 (m’), ifpsr. 
SP) (m) = 0 (mB), if B>r. 


Also, 
* D,” [a, (x)] = 0 (x-#*4), 


for any integer r > 0, as x > 9, 
So, if y > 7, 
S™ (m) - ay, (mt) O (m- + 47) 
dmy 


-0 (1), if > y— 5. 





*Su ()=0(T) as xX > co, 


1 j J _ 1 F 
So ap (x)=O (ari): Da’ @u@)= 2 PO) De! In + De? S [<<]: But D,* Jp can be 
expressed as a linear combination of Jy’s with the index » ranging from u—s tow+s 


Therefore, Dx Ju=O (T;) 


Hence Daf Jn» Da? * [=] va (arpe=)= ” (sxz3) 
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If y <!, 


(y) 
S™ (m) a a, (mt)= O (m- + n+7) 


=o(l1), ifu>r—4. 
Hence, from (3-4), we have 


$ a, a, (nt)=0(1)+(— 1+? f s+) (x) D,4+ % a, (xt) dx, 
0 


if wp >h+h. (3-5) 
Now, (3) ; 


St? +) (m) = S+ FFT=P) (m) = wes 4 af (m— x)"~* S$ (x) dx, 


so that 


oO 


7 ] ; hs oe , es 
-raeEca Jf D a, (xt) dx fe u) SOX u) du. 
° 0 


: 1 " (7) - ee h a f h +2 o 3 
P(U+h—n) f S™ (u) du f (x— u) D. Oy (xt) dx. (3 6) 


u 


From (3-5), (3-6), we therefore obtain, 


2 dy Oy (nt) 
0 


00+ oY, se (u) tu Fewer D,!*2 a, (xt) dx, (3-7) 
provided 


m 


f S (u) du f (x— uy" otet D,’ +2 a, (xt) dx >O0asm—+co. (3-8) 
0 m 
Now, for a fixed ¢, we have by the Second Mean Value Theorem, 


| fe- uf’-* DP +2 a, (xt) dx 


< (m—u)"-” Max 


m > m 


-(m— u)’-7 0 (m-#- 4) (3-9) 


f D,’+? a, (xt) dx 


m 
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and 


=m F | 80 (u)| Gm wy =" d= O (ol 2-9 (3-10) 


oO 


=0(1), ifu>h +4. 
From (3-9) and (3-10), we see that (3-8) is satisfied, and so, 


5 a, 4, (nt) 


= (”) a oe. D,’ +? t) dx, 
0+ po 5S wo au fc u) a, (xt) dx 


if p>h+h. 


Letting m — oo, we have, 


s 
2 a, 4, (nt) = 
0 


co 


(— 1 *  Panas ee (x—uy’-”- D+ 2a, (xt) dx, (3-11) 
r(1+h—r) r ‘ | 


provided the right side exists. 


In order to prove the theorem, we have only to prove that, if 


co 


co 
@ (t) = f S™ (u) du fie- af-*- 4+ * a, (xt) dx, (3-12) 


then 

(i) ® (1) is finite for fixed ¢ > 0, (3-13) 
and 

(ii) & (t) >Oast—0. (3-14) 


Setting x= uy in (3-12), we have 


co 


(1) =1'+ bd (u) dur uw we (ae 1)’-*%+ a+ 2) (yut) dy, (3-15) 


v 


where, 





ou Out) = aemay [xt Ow) | 











0) 


x, 


14) 


15) 
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Now, let 


I(z) = f (y— 1-7 a+ % (yz) dy 
» | 


so that, from (3-15), we have 


co 


@(1) = t+? f $7) (u) dus w+ *=" (ut) du 
oO 
We shall evaluate the order of I (z) 
1+ 1/z co 
I (z) =f b f = IJ,+ I, say. 
5 1+4/z 


For |z| 21, 
1+i/z 


lhis f [\(y— 1)*-” a,** 2 (yz) | dy 


1 
1+1fe 


< O-*-h f (y—1)'-" dy 


1 


= O(z- +") p+ r), 
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(3-16) 


(3-17) 


(3-18) 


(3-19) 


(3-20) 


And, 
co 
ire |(y— 1*-r a4 + (yz) | dy 
1+1/z m 
™ | } , 
E a,(* +) (yZ) (y- 1-7] 
1+ 1i/s 
co 
nies sae. a, +1) (yz) (y—- 1)*-*-1 dy 
= O (z~(#+b-4+7) 1+ 1/z 
For |z| < 1, 


1+ 1/2 


[il = f (y—1)*-* a,4*% (yz) dy 


= 0 (z-4-7+0) 
[I,] =O @--7+ 9), 
Hence, from (3-18), (3-19), (3-20) we have, for |z| = 1, 
[1 @)|= 0 G-M+tb-A07 
and, for |z| <1, from (3-18), (3-21), (3-22), 
[I (z)|=O (2-@-7+ 9), 


(3-21) 
(3-22) 


(3-23) 


(3-24) 
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Now, from (3-17), we have 
1 


/t oe 
@ (t) =t* +2 [ f S”) (uv) wt+4-*% T (ut) du+ f] 


1/t 
_ K, + K,, Say. (3-25) 
1/t 
| : 1 
a h + 2 1 hk eee oe < a oy 
K.= 1 f O(l)u O (Gayrrei) since 0 <U< 5 (3-26) 
° 
K, {? 2 f 0 (i) ul A O (cera) du, since ut > Lt. (3-27) 
P) 


Now, for a fixed ¢, K, is finite, and K, exists if 1» > r+ 4, so that ® (1) exists 
for a fixed ¢. 

It is easily seen on the other hand that, if we let t—0, K,-0()), 
K,= 0 (i), if only » > r+ 4, so that 6(4) 0 as ¢t > 0. 


Hence, from (3-11), we observe that 
co , 
b, (1) a) Gy ty (nt) 
n=0 
converges for ¢ > 0, and ¢, (1) > 0 as t> 0. 


Therefore 2 a, is summable J, for w > r+ 4. 


Remarks. As a generalization of the above theorem we can prove that 
if a series = ay is summable (C, a), then it is summable (J * n) to the same 
sum, provided K(u+4)>a+1, where K is an integer greater than zero. 
The proof follows on exactly similar lines. This theorem includes, as a 
special case, the following theorem of Kogbetliantz*: 


If = a, is summable (C, a) to the sum S, then for p > a+ 1, we have, 


for, we have only to set K = p, »= 4 in the above theorems. 


Theorem I is ‘ the best possible result’ as can be seen from an example 
in Fourier series given by Paley.’ 

§4. THeorREM Ii. [f2 a, is summable J a» and (i) 2 a, oy (nt) converges 
in such a way as to permit integration term-by-term, for every t> 0, 

(ii) Za, a, (nt) = d, (t)= O(1), as t-> 00, then, Z a, is summable (C, r) 
forr>pt+h. 
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Proof. Without loss of generality the J, sum may be taken to be zero 


co 


Let F (eo)= o* * fa, (ot) dy ( 0 1 dt, (8-1) 
0 
Then we shall first prove that 


Lim F (w)= Lim 4, (f) = 0, provided v > 
w> co t-0 


> let 


bdo) YW 


Given « > 0, | dy (| < ie for |t| < 6(«), where K is a constant to be 
defined presently. 


Now 


6 (e) 


F (w) oi | 


fof 


6 (e) 
I,+ I., say. 
5 (€) 
[Ll <p: o™ a4 | a, (wt) | 242 de. 
cu 
w96 (t) 


K x*#+1 1a, (x)| dx. 


co 
€ 


<x f x1 Ja, (x) | dx 


(8-3) 
° 
< «, if y> M+ 14, 


co 


where fx t |} a, (x)| dx = K, since a, (x)= O (x-%*#) 


0 


co 


\I,| < 7 | | a, (x) | x +2 dx 


o 


< §; 


where 


|b, (0) | < A, for, given 3 we can choose w so large that 
(8-4) 
dw 
Hence, from (8-2), (8°3), (8°4), we have 


F (w) =0(I) (8-5) 
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co 


F (w) = 2 a, w+ + fm “1 a, (wt) a, (nt) dt (8-51) 


0 
co 


wt 2 f p2m+t a, (wt) a, (nt) dt 


° 
co 


; eet J, (wt) J, (at) 
— 9 tv = = 2 € Ma a es cies 
2+ Put 1+) f (cot? (nt dt 
a tr CF, (ot) Ip (at) 
nP jp-P -1 
0 
where K is some constant. 


-K- 





From the theory of Bessel functions,® we have, 


co 


f= 3,(bt) b I (du+ 4v—4A+ 4) 





tr ~ W@-*tI T+) PGA+ Ww — rth) 
+v—A+1 v—A-p+l 
ee es 


_ 


0 


provided,0O<b<a, R(u+v+1)>R(A)> —-1. 
From (8-7) we deduce that 


co 


f J, (ot) In Ot) . (1- _— 
WwW 


ge-e-1 Qv—-K-1 wth+2-¥ I'(v—p) 





° e 
if0<n< aw, 


and =0, ifn > w. (8-8) 
From (8-5), (8:5 1), (8:6) and (8-8), we observe that 
\y-K-B-1 
K-F(w)= 2 a, (i- al -> Oas w-> co (8-9) 
a) wo" 
where K is a constant. 

Hence 2 a, is summable (R; n?, v—p— 1) for v>2n+ 14. By the 
second theorem of consistency for summable series due to Hardy,? 2a, 
is summable (R; 7, v—y» +1); that is, 2a, ’is summable (C, r), for 
r>pt t. 

Remarks. As a generalization of Theorem II, we can prove that: 


If (i) the series 2 a,, a, (A, t) converges for all t > 0, and uniformly 
in every finite interval, and 





Bessel Summation of Series 
(ii) , (¢) is bounded as ¢t > co, and 
(ili) ¢, (t) > 5+ co as t> 0, 
then 2 a, is summable (R; A, K), for K > p+ $. 


It must be stated, in passing, that Theorem II can be put in the 
following form in which it finds application usually, viz., 


If X a, is summable J,, then it is summable (C, r) for r > w+ 14. 


This follows easily, since a,—0(n*+#), and so Ya, a, (nt) converges 
uniformly and absolutely for »>p+ 1. But thisis not the best possible 
result. The restriction on r seems to admit further relaxation. 


However, Theorem II is ‘the best possible result’ as can be seen from 
the classical example of Hahn? in Fourier series. 
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